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' Inspired by the need for effective stochastic models to describe the complex behavior of 

biological motor proteins that move on linear tracks, exact results are derived for the velocity 
and dispersion of simple linear sequential models (or one-dimensional random walks) with general 
waiting-time distributions. The concept of "mechanicity" is introduced in order to conveniently 
quantify departures from simple "chemical," kinetic rate processes, and its significance is briefly 
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indicated. The results are extended to more elaborate models that have finite side-branches and 
include death processes (to represent the detachment of a motor from the track). 



I. INTRODUCTION AND SUMMARY 



Motor proteins such as kinesins, dyneins, myosins, DNA and RNA polymerases, are important for the biological 
functioning of cells. Consuming energy obtained from the hydrolysis of ATP or related compounds, and moving along 
rigid linear tracks (microtubules, actin filaments, DNA molecules, etc.), they play significant roles in cell division, 
cellular transport, muscle contraction and genetic transcription 1 ' 2 . Such molecular motors can move with velocities, 
V, up to 1000 nm/s 3-6 and may sustain an external load, F, of 5-8 pN for kinesins 6 and up to 30-40 pN for DNA 
and RNA polymerases 7 ' 8 . Understanding the detailed mechanism of the functioning of motor proteins is a major 
T [ challenge of modern biology. 

In recent years, significant advances have been made in experimental techniques for studying motor proteins: one can 
now observe and investigate accurately the mechanical properties of single molecules over wide parameter ranges 3-8 . 
However, our theoretical understanding of how these proteins work is still incomplete. Theoretical modeling of 
molecular motors has followed two main directions. One class of models is based on "thermal ratchets" 9-12 in which 
the motor protein is viewed as a Brownian particle that diffuses in two or more periodic but asymmetric potentials 
between which it switches stochastically. Another, more traditional "chemical" approach is based on multistate kinetic 
descriptions of the motion with various rate processes determining the transitions between the states 13-23 . In this 
paper we consider various extensions of the chemical kinetic schemes for which we derive exact, closed form results in 
>^ | terms of the underlying transition rate parameters. 

. The simplest periodic sequential kinetic model assumes that a motor protein molecule moves along a periodic 
molecular track (a microtubule in the case of a kinesin) and binds at sites at x = Id (I = 0, ±1, ±2, • • •). It is 
postulated that there are N discrete states, j = 0, 1, • • • , N — 1, between consecutive binding sites and the motor 
protein in state ji ("at" site /) can move forward to state (j + l)i (normally pulling a load) at a rate Uj and can slip 



backward, to state (j — 1)/ at a rate Wj 18 - 19 . This basic iV-state model is, clearly, precisely isomorphic to a discrete 



3 

but, in general, biased random walk on a periodic one-dimensional lattice. Such random walks are of broad interest 
for a variety of applications 24-29 , in particular, for studying diffusion in random environments 28-30 which may, for 
• example, be approached by allowing the period N to become infinite 28 . Indeed, some time ago Derrida |28| presented 
a mathematical approach that provides formally exact and explicit formulas for the asymptotic drift velocity 

o' 

V = Vodu^wj}) = Km ^(x(t)), (I) 



dt 

and for the dispersion (or effective diffusion constant) 



5_i 

D = D ({ Uj , Wj }) = \ Urn j t [(x 2 (t)) (.*(i)} 2 ] , (2) 

where x(t) is the spatial position of the motor protein molecule or, equally, of the random walker along the linear 
track at time t. Derrida's results have recently been exploited in developing theory for molecular motors 17-20,23 . A 
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valuable feature of the exact closed expressions is that it is quick and straightforward to explore the precise effects of 
a wide range of parameter values, revealing a variety of different characteristic types of behavior 18-20 , in contrast to 
employing approximate numerical schemes or Monte Carlo simulations, or to studying over-simplified models . 

The "purpose" of a motor protein in a biological setting is to move various biochemical entities, such as vesicles, etc., 
that exert load forces, F, which act on the protein molecule. The questions of what force can be exerted by a motor 
protein within the basic model and how the transition rates should be changed by the load, F, have been discussed 
critically 18,19 . The resulting calculations for these models show qualitative and good semi-quantitative agreement 
with the experimentally observed behavior of normal two-headed kinesins 3-8 ' 18 ' 19 ' 23 . 

In order to describe motor proteins more realistically, however, and to understand the behavior of proteins other 
than two-headed kinesins 7 ' 8 ' 31 ' 32 , extensions of the basic A-state periodic sequential kinetic model were recently 
introduced 20 . These extended models take into consideration the complexity of real biochemical cycles (multiple 
kinetic paths and branched states) and allow for the irreversible stochastic detachment of molecular motors from 
the linear track which is always observed. Exact expressions for velocities and dispersions have been derived for 
these extended chemical kinetic models 20 in terms of the basic forward and backward rate constants, Uj and Wj, the 
corresponding rate constants for the side branches, and the death rates from any of the N intermediate states. 

Now, the basic concept underlying chemical kinetic models is the idea that the motion of the motor protein is 
essentially "chemical," i.e., the molecule undergoes a transition or "jump" from one chemical state to a nearest 
neighbor state (in the kinetic diagram) at a given rate having "forgotten" how it arrived at the state. The time 
intervals between such jumps are thus distributed exponentially according to Poisson statistics. In other words, the 
time intervals are described by exponential waiting-time distribution functions in which the coefficient in the exponent 
represents the overall rate of transition from the state in question. 

A conceptually possible alternative picture of the motion is that it is "purely mechanical" so that the motor 
works like a clock, jumping from state to state within narrowly distributed time intervals. Real molecular motors 
work with a surprisingly high mechanical efficiency (which may be estimated as 80-90% for kinesins). This, in turn, 
suggests that some of the steps in the dynamical sequence of the motion might reasonably be described as rather more 
mechanical than purely random in nature. In that case, the waiting-time distribution functions for the time intervals 
marking transitions between different internal states might be modeled more effectively or more economically 23 as 
non-exponential. To that end we present here a generalization of the basic periodic chemical kinetic models, and of 
their extensions, that explicitly considers general waiting-time distributions. By using the results of Montroll and 
coworkers 24-27 we are, again, able to obtain explicit, exact results for V and D for the various models. We also 
explore briefly the quantitative significance of the "mechanical" factors that could represent departures from purely 
chemical kinetic processes in the motion of molecular motors. 

To proceed, consider first the extension of the simple sequential kinetic model (with period N) shown schemati- 
cally in Fig. 1. To characterize the dynamics of the motor protein (or random walkers) we introduce waiting-time 
distribution functions ipj(t) and ipj(t), where ipj~(t)dt is the probability of jumping one step forward from state j 
in the time interval t to t + dt after arriving in state j, while ipj (t)dt is the corresponding probability of jumping 
one step backward from state j, and ipj(t)dt is the probability of attempting to move in the same time interval but 
failing to do so and hence restarting at the same site after the attempt. The allowance for tp®{t) is potentially useful 
in modeling motor proteins since it can describe, say, the "futile" hydrolysis of an ATP molecule without achieving a 
forward step. 

A more complicated model allows for the possibility of irreversible detachments or "deaths" . Deaths are described 
by a distribution Vf (*) sucn that ipj(t)dt is the probability of leaving from state j in the time interval t to t + dt but 
not appearing in any other state 20 . Alternatively, the detachment or death process may carry the motor, or walker, 
to a reservoir (or graveyard) from which no returns are allowed. 

Finally, another possible extension of the simple kinetic sequential scheme is to incorporate finite side-branch 
processes 20 : see Fig. 2. If state (j,k)i labels a state A; on a branch emanating from the primary state ji, the 
dynamics can be described by additional distributions ipjk(t) where ipj k (t)dt is the probability of jumping one step 
further out from the branch state (j, k)i in the time interval t to t + dt (k = 0, 1, • • • , L — 1), while ijij k {t)dt is the 
corresponding probability of jumping one step back in towards the primary state ji = (j, 0); from the branch state 
(j, k)i (k = 1, 2, • • • , L). Note that without loss of generality we may assume that all side branches are of the same 
length L: see Ref. 20. 

The waiting-time distribution functions respect the basic TV-state periodicity so that 

4 k (t)=^ j±Ntk (t), (3) 
where ^ = +, — , 0, S, /?, or 7. Furthermore, for any state (j, k) normalization requires 
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(4) 



For all these extensions of the original chemical kinetic model we have derived explicit general expressions for the 
drift velocities and for the dispersions by generalizing Dcrrida's method 28 . When appropriate exponential waiting- 
time distribution functions are substituted, we recover Derrida's original formulas 28 for the simple linear model and, 
likewise, our previous results for the extended kinetic models 20 . For convenience, we present all our concrete results, 
namely, the expressions for the mean velocities and for the dispersions, in this section. The detailed, and unavoidably 
somewhat involved, calculations and derivations are described in Sections III-V. 

To quantify the importance of "purely mechanical" factors in the motion of motor proteins the concept of mechanic- 
ity is introduced and discussed in Section II. The mechanicity varies from for a "chemical" or Poisson process (with 
exponential waiting-time distribution functions) to 1 for a purely mechanical process (with clock-work or delta-function 
waiting-time distributions). It provides a convenient quantitative measure of deviations from a simple chemical picture 
of the dynamics of molecular motors which may well prove useful in applications 23 . 

To simplify the presentation of our results it is convenient to introduce the overall waiting-time distribution, ipj,k(t), 
which is merely the sum of all the distinct distributions associated with the state (j, k). (Recall that k — corresponds 
to primary states, ji = (j, 0);, on the linear sequence). We will need the Laplace transform of the overall distribution 
defined by 

poo 

i>j,k( s ) = / e~ st iij, k (t)dt 
Jo 

= 1 - S/J.x-j.k + \s 2 ll 2 : 3 .k , (5) 

where ^i n -j.k (n — 0, 1, • • •) is the n-th moment of iftj t k(t) and we have used the normalization condition (4) to conclude 
fJ<0;j,k = 1- The transforms, ipj k (s), of the partial waiting-time distributions and their moments jjL n .^ k , are defined 
similarly. Of course, for the specific models that we consider some of the ipj k will be absent (and, so, can be set to 
zero). 

The analysis reveals that a crucial role is played by the associated waiting-time rate distributions, ipj k (t), which 
are defined via their Laplace transforms 

= - (6) 

where (, as before, can be +,—,0,5,0 or 7. Note that it is the transform Vj,fe 01 the total or overall waiting-time 
distribution that appears here in the denominator. 

It then transpires, as seen below and proved in Sections III-V, that in terms of the <fj k (s) one can readily define 
effective kinetic transition rates: Uj and Wj , for the forward and backward rates along the primary linear sequence; 
Sj for the death processes from states j = (j, 0); and side-branch rates, (3j t k and jj ik , outward and inward from the 
states (j, k). In the main, especially for the velocity V^[{?/;j fc (t)}], these effective kinetic rates play the same role in our 
results as do the simple (Poissonian) rates, denoted by the same symbols 18-20 , in the analysis of the original kinetic 
models with exponential waiting-time distributions. Furthermore, for the velocity and dispersion (in the absence of 
death processes), only low order moments of the ipj k (t) enter explicitly. These, in turn can be expressed in terms of 
low order waiting-time moments, fi n -,j,k and ir n .- k . To be explicit, dropping the state labels (j, k), we can write 

fi(s) = v% - svi + \s 2 vi , (7) 

where the rate moments u£ are expressed in terms of the \x n and /i£ in Table I. We will assume that all the moments 
displayed in the Table are finite. 

As found previously in analyzing the various simple kinetic models, our expressions for V and D in the extended 
models with waiting-time distributions depend on certain linear sequential products of rate ratios. Following 20 , we 
thus define two types of product for the sequential model with waiting times, namely, 

n^n^i and nf = n^±l = ^±inf. (8) 

J J - J - Ui J J -- L Ui Wj J 

The model with branches requires one more type, namely, the branch products 
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k 

ii a frfcl. (9) 

3 11 /V' ' 

2 — 1 J ' 

For the model with deaths or irreversible detachments we introduce modified analogs of (8), namely, II^ and IX- fe , 
that are obtained simply by substituting Uj and Wj by the "renormalized" values 

u j =u j e j+1 /e j and w, = UjEj-i/ e h (10) 

where the periodic renormalization coefficients, Ej = Ej±N {j = 0, 1, • • • , N — 1), arc conveniently normalized by the 
condition 

(11) 

The remaining coefficients Sj then form the components of the right eigenvector, e = [ej], corresponding to the 
smallest eigenvalue, A = \({uj,Wj}), of the NxN transition rate matrix M[{uj, Wj; 5j}] which is defined by the 
nonzero elements 

M i,3-i = ~Wj, M jtj = Uj + wj + Sj, M jJ+1 = - Uj , (12) 

and, because of periodicity, 

M ,_i = M ,jv_i = -w and Mjv-i.jv = Mw-i,o = ~Un-i- (13) 
With these preliminaries established we can now present our explicit results in fairly compact form. 

A. Results for the Simple Sequential Model with Waiting Times 

The formal expression for the velocity in this case is precisely the same as derived by Derrida 28 for the sequential 
chemical kinetic model, that is 

V = d(l - Ti?)/R N , (14) 

where, using the definitions given above, 



N-l 

3=0 



N-l 

-tj+k 

fe=l 



(15) 



and d is a distance between neighboring binding sites on the linear track. However, the effective transition rates Uj 
and Wj are now given by 

~+ + JT^ (*)<** 

Uj, Ws= v. {8 = 0) = Vo; . = (16) 

where we have used (7) and Table I and may recall that ipj = xp^ + tpj + ip®. 
The dispersion takes a new, more complex form which can be written 

D = D + D 1 , (17) 

where the first term is given by 

Do = (d/N){[VS N + dU N ]/(R N ) 2 -±(N + 2)V}, (18) 

N-l N-l N-l 

S N = S 3 E ( fc + fyk+j+U U N = ^ U 3 r 3 S 3> ( 19 ) 

j=0 k=0 3=0 

in which new coefficients Sj are determined by 
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Sj = Uj 1 



JV-1 

i + E njif 



fe=i 

The second contribution to the dispersion is found to be 

JV-l 



(20) 



D 1 = (d/N)[NV/(R N f\ J2 *M r i 97+irj+i), (21) 

3=0 



where new, "nonexponential parameters", gf, are defined by 

/>OC 

/ t<pf(t)dt = -vf. j . (22) 
Jo 



dtpf 



ds 



s=0 



For the relation of these parameters to the waiting-time moments, see Table I. When the waiting times are expo- 
nentially distributed with, in particular, ipj~{t) , ipj~{t), and ipj(t) each proportional to exp(— Cjt), the parameters gf 
vanish identically. When this occurs for all j, the dispersion D is given by Do alone and, as mentioned, Derrida's 
original formula is recovered 18,19 ' 28 . We remark that the g^ are typically negative but may, in fact, also be positive: 
their character is discussed in more concrete terms in Section II. 

B. Sequential Model with Branches and Waiting-Times 

For the models with branches of finite length (see Fig. 2) the velocity is given by 

V/ /3 (/3, 7 ) = rf(l-nf)/<, (23) 
which is identical to the expression (14) for the unbranched models except for the modified functions 

JV-l L 

<=E^< rf=r J [l + j:nf*] > (24) 

3=0 k=l 

in which the effective transition rates in and out of the branch states are defined in parallel to (16) by 

= = 0), 7i,k = £j, fc (* = 0). (25) 

It is appropriate to recall however, that the rates fUjfi, u 3 - and Wj involve ipj,o(t) which now entails the three distributions 
■0+, and tp? . 

The dispersion can be expressed as 

Dp = D 0>/ 9 + Di,/j + D 2t0 , (26) 
where the first two terms are very similar to D and D\ for the unbranched models: explicitly we find 

D ,p = (d/N){[VpS N + c<]/(<) 2 -\{N + 2)V P }, (27) 

N-l 

D U3 = (d/N)[NV /(R N f} E s M r i 9]+ir j+ i), (28) 
where the modified functions, analogous to Sn and Un in (19), are 

JV-l JV-l JV-l 

s n = E s " E ( fc + ^Vi' u n = E w?, ( 29 ) 

j=0 k=0 j=0 

in which, in parallel to (20), we require 
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JV-1 



z=l fc=l \ z=l / 

The third term in (26), which is a contribution due to the presence of branches, is given by 

JV-1 L 



where, for the coefficients, we have 



3,k 



0=0 fe=l 



2^ J 9 3,k-l LL j + 9j,k LL j 



fc-1 



+ E§7 n f ( E nf'-^.nf^ + ^n; 
i=i V*=fc-i 

The nonexponential parameters, fc and <?J fc , are defined in precise analogy to (22) by 

9j,k = {d^l k /ds)\ s=0 = -vl tjjk , g] k = (d^] k /ds)\ s=0 = -vj. jjk , 



0,k-i 



(30) 



(31) 



(32) 



(33) 



where, again, the relation to the waiting-time moments follows from Table I. Similarly, if ipj k (t) = k cxp(— Cj^t) 
holds for all (j, k), the results of Ref. 20 arc once more obtained. Section II discusses expressions for the g^ k in terms 
of associated "mechanicities" Mf . . 



C. Sequential Model with Deaths and Waiting-Times 

The mean velocity is now given by 



Vs = d(l - Tl?)/R 



JV> 



(34) 



where II^ and the fj (appearing below) are defined, using Uj and Wj [see (10)], in exact analogy to (8) and (15). 
However, in contrast to (16), the effective transition rates Uj and Wj and the death rate Sj, which enter the expression 
for R S N , require new definitions. Specifically, we find 



<p?(a = -A) 



/ °°(e+>*-l)^(t)dt' 



(35) 



and 5j = !pj(s=— A) which can, likewise, be expressed in the integral form exhibited in (35), where it may be recalled 
that ipj = ipj' + tpj + ipj + tpj . In these expressions A is, as stated above, the smallest eigenvalue of the transition 
rate matrix M[{u,-, Wj; Sj}] [see (12) and (13)]. The function R S N is then given by 



R 



N 



JV-1 
3=0 



r 5 j = [l + (l-a j )g 5 j }r j , 



where the modified death and non-exponential parameters g$ and gf are embodied in 

a 3 = [gt( s 3+T- - £ j)+9j~(£j-i - £ 3 )]/( £ 39j), 



while the parameters themselves are now given by 



„<5 _ 



(dtf/ds)\ s 



-A > 



(d$f/ds)U 



(36) 



(37) 



(38) 



where the integral expressions corresponding to that in (35) are somewhat more elaborate but may be found strait- 
forwardly from (6) and will mirror the form of v£ in Table I. Recall from (10)-(13) that the coefficients Sj form the 
right eigenvector of M for A, the smallest eigenvalue. 
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Finally, the dispersion can be written as 

D 5 = D , S + D h s + D 2 ,s, (39) 

where the first term is given by 

Do.s = (d/N){[V s S 5 N + dU s N ]/(R s N ) 2 -±(N + 2)V S }, (40) 
in which Sfj and U${ are defined in precise analogy to (19), that is, 

JV-l N-l N-l 

s n = T, s j E ( fc + i y k+]+ i> ^ = E wl> ( 41 ) 

j=0 fc=0 j=0 

fj being defined just after (34), while the analog of (20) is 

N-l 



1 + (1 - a-)g\ + £ [1 + (1 - a j - k )g s j _ k ]IL)i- 1 k 
fe=i 

The second term in (39) is then given by 



(42) 



JV-1 

D liS = (d/N)[NV s /(R s N f] ]T * S M?i gj+irj+i), (43) 

j=o 

where we have introduced the renormalized modified nonexponential parameters 

gf =gfs j ±i/s j . (44) 
Finally, the last term in the dispersion, which arises solely because of the possibility of death or detachment, is given 

by 

N-l 

D 2 , s = \{Vi/R%) £(1 - Mhfc, (45) 

1=0 

where the new functions 

Pi = H(e j+ i - ej) + hj (£,_! - e^/iejh*), (46) 
incorporate the second-order nonexponential parameters 



roc 

/ t 2 e xt ^f(t)dt, (47) 
Jo 



hf = (d 2 lpf/ds 2 )\ s= . x = / t 2 e Xt <pf{t)dt, 

Jo 

/>OC 

h 5 j = (d 2 $/ds 2 )\ s =- x = / t 2 e xt ^{t)dt. (48) 

Jo 

It is interesting that these higher moments of the waiting-time rates, ^pf{i) and <Pj(t), arise only when death processes 
come into play. They can be written as integrals of the corresponding ipj(t), as in (35), but will then exhibit the 
structure of i>2 m Table I and entail, in particular, the modified third moments / °° t 3 e xt tp^dt, etc. Note that when 
the probability of death or detachment vanishes, one has A = and the coefficients Uj , (3j ,h^,h s ^ and g 5 - all vanish; 
then we recapture the results stated above for simple sequential model with waiting times. 
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II. DEGREES OF MECHANICITY 



Our analysis of the velocity and dispersion for the stochastic models with waiting-time distributions has revealed 
that in terms of the effective rates Uj, Wj, etc. [as defined in (16), (25), and (35)] deviations from Poisson processes 
with nonexponetial waiting times do not change the formal expressions for the velocities: compare, e.g., (14), (23) 
and (34) with the corresponding results presented in Ref. 20. On the other hand, the expressions for the dispersions 
change dramatically and, in particular, involve new "nonexponetial parameters" like g^~ [see (22), (33), (39), etc.] 
that evidently provide a measure of some sort for the departures from simple "chemical," kinetic processes. In order 
to gain a more concrete and intuitive picture of what these departures involve, we introduce a quantitative concept 
which we call the mechanicity, Mj, of the specific transition process described by the waiting-time distribution ipj(t) 
(where, as previously, ^ = +, — , 0, • • •) . 

The basic idea is to discriminate, in an explicit way, a standard exponential distribution from distributions that 
depart from it in varying degrees, and from a sharp, "clockwork" distribution of zero width. Now if t is the mean 

— — 2 

waiting time for a particular process, the mean square deviation (At) 2 = t 2 — t provides a natural measure of the 

— 2 

width of the distribution. Then the dimcnsionlcss ratio O = (At) 2 /t represents a scale-free index of the relative 
width. In the purely mechanical or clockw ork limit At = t — t must vanish identically and so O = 0; conversely, if 
the distribution %j)(t) yielding t and (At) 2 is a simple exponential, one has = 1. Thus the mechanicity parameter 
M = 1 — O vanishes for a Poisson or "chemical" process but attains the value unity for a purely mechanical process. 
More illustratively 18 , suppose the waiting-time distribution has the familiar general form 

${t) = Qt v - X e~ at [v > 0). (49) 

An elementary calculation then yields 

M = 1-6 = (50) 

By construction, the chemical limit is described by v — 1; conversely, the mechanical limit is realized when v — > oo. 
This example also shows that M may be negative (and, then, indefinitely large). Evidently, this arises when ip(t) is 
sharply peaked at the origin exhibiting a power law behavior ~ l/t x with \ = |M|/(1 + \M\) (for M < 0). 
In applying these considerations to a general waiting-time distribution, ipj(t), with moments 

poo 

<*C>i=*4-= / *"^-(*)*. (51) 

JO 

a little care is needed since the zeroth moment, fj^.j = (t®)j, represents the total probability that the transition of 
type £ occurs which, in general, is less than unity. The appropriate definition of the mechanicity for the process £ 
from a state j is thus 

Mf = l-Q< = 2-{t° c ) j (1? c ) j /(t c ) 2 j . (52) 

Now in applying the various expressions displayed in Section I for the dispersions, D, it would be convenient if the 
nonexponential parameters, g^ (or g? and gj), could be expressed directly in terms of the associated mechanicities, 
Mf (and/or and Mj) and the effective rates 

U 3 = (t+)j/ T j and W 3 = (t-)j/Tj (53) 

(and/or fi^k and 7j,fc), where if, for simplicity, wc suppose only forward and backward processes act from state j, the 
mean dwell time is just 

Tj = (t+ )j + (t-)j = (uj +W,)- 1 . (54) 

However, because the g^ k are defined via the rate distributions V'jfc(^) which in turn, as seen in (6), require the 
total waiting-time distributions, tpj^it), matters are not entirely straightforward. If one stays with the simplest case 
= ip? = ipj = 0) a single extra parameter proves essential: this might, for example, be taken as the dimensionless 
ratio 

^ = (Uh/Tj < 1. (55) 
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But the resulting formula, derived from 



at = [(tlh + (t-h W+hfa] . (56) 

which in turn follows with the aid of Table I, has a paradoxical feature (for which reason we do not quote it). Namely, 
even if the separate mechanicities, Mf and M~ , vanish the parameters gf and gj do not, in general, vanish! The 
reason is that in a pure kinetic scheme with forward and backward rates Uj and Wj describing departures from the 
same state j, the resulting waiting-time distributions, tl>f{t) and share a common exponential factor, namely, 

exp[— (uj + Wj)t] = e^ t / Tj . Thus even if ipj~(t) an d V'jW are both simple exponentials, the overall process will not 
have a simple chemical description unless the + and — decay rates match. More generally, however, if Mf = MJ = 
and the (single) additional condition 

<4>iH*±W(<t + >? + <*->?) (57) 

is met then, indeed, gf and gj vanish. (Similar considerations apply, of course, to the behavior on branches; but the 
extra process available at a primary state, where a branch starts, must not be forgotten.) 

Despite these conceptual complications, one can devise instructive examples with simpler dependence on the me- 
chanicities. One useful case when only forward and backward transitions from a state j occur is described by 

^(*) =Qft v± ~ 1 e~ v±t/rs , (58) 

where, with (53), the relations 

Qf = ( Uj , Wj )r]- V± v^/T{u ± ) (59) 
ensure that the normalization condition /q°°(V'/ + ipj)dt = 1 is satisfied. The mechanicities are clearly 

M+ = 1 - vZ 1 and MJ = 1 - vZ 1 , (60) 
where, of course, v + and v_ could also depend on the state j; the nonexponential parameters are then given by 

g f = _i( Uj) w 3 )rf( Uj M+ + Wj M7). (61) 

This result does depend only on the effective rates and the associated mechanicities and the gf do vanish when 
— > in accord with the naive expectations. 
Similar examples can be devised when branching and death processes occur from a state j. For modeling purposes 23 
we believe the associated loss of full generality is likely to be insignificant. But note that, by (61), the gj are always 
negative for the special cases considered; if, however, v + = v- = 1 so that M+ = MJ = but there are distinct time 
constants, r+ ^ tJ , then gf and gj will have opposite signs. We stress, nonetheless, that the results presented in 
Section I apply for quite general waiting-time distributions: the mechanicity may be regarded as an auxiliary concept 
of intuitive and descriptive value. 



III. PERIODIC SEQUENTIAL MODEL WITH GENERAL WAITING TIMES 

To derive the results presented in Section I, consider first the general periodic sequential model with waiting- 
time distributions as specified in Fig. 1. This model can be regarded as a one-dimensional continuous-time random 
walk with N internal states (j = 0, 1, • • • , N — 1), a class of walks considered some decades ago by Montroll and 
coworkers 24 ^ 27 ' 29 . The crucial result, demonstrated by Landman, Montroll and Shlcsingcr in 1977 27 , is that the 
probability Pj(l,t) of finding the walker at site I in state j at time t satisfies the generalized master equation 

j f Pj (I, t) = J o { <pU (T)Pj-i (I, t-r) + ^J +1 (r)P j+1 (I, t - t) 

-[<p+(T)+ i pJ(r)]P j (l,t-T)}dT, (62) 

where the relaxation or memory functions, <pf(t), are related directly to the waiting-time distribution functions, 
ipf(t), precisely as specified in Eq. (6). This master equation replaces the simple kinetic rate equations which were 
the starting points of the previous analyses 17-20 ' 28 . 
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For our purposes, we may, without loss of generality, assume that the initial condition is Pj(l, 0) = P^Si^: i.e., the 
walker starts at the origin x = I = 0. Conservation of probability then dictates 



+00 JV-1 

]T £P,-(i,t) = l (aUt). 



(63) 



On the other hand, in any arbitrary state ji in the kinetic picture the normalization requirement (4) yields 

/ ty+(t)+^7 (*)+$(*)]<** = !, (64) 
Jo 

or, in terms of Laplace transforms, 

$f(8 = 0) + $j(s = 0) + ^(s = 0) = 1. (65) 

In order to find the drift velocity V and dispersion D we now generalize Derrida's method 28 by defining two auxiliary 
functions for each state j, namely, 



(66) 



(67) 



B,-(t)= J] W), C o(t)^ E (3 + Nl)Pj(ht)- 

l — — oo l — — oo 

The generalized master equation (62) then yields 

jB.it) = J {p+^B^it -t) + <pJ +1 { T )B 3+1 {t - r) 
-[ l p+(r)+ i pJ(r)}B j (t-r)}dT. 

Similarly, we obtain 

jC 3 (t) = J* {^-i (r)Q_!(t - r) + <pJ +1 (T)C j+1 {t r) [pf(r) + <pj{r)] C 3 (t r) 
-<pJ +1 (T)B j+1 (t - r) + iff^Bj.^t - r)} dr. 
Again following Derrida 28 we introduce the ansatz 

which should be valid when t — > oo. Because of the periodicity in j we have 

bj+N — bj, a j+N — and Tj+jv = Tj. 
After long times a steady state, dBj/dt = 0, will be achieved. Then, recalling in (35) that J °° ip^(t)dt = <pf{s = 0), 



(68) 



(69) 



(70) 



we introduce the effective transition rates Uj and Wj dchncd, in anticipation, in (16). Thus (67) yields 

= Uj-ibj-i + Wj + ib 0+ i - (uj + wj)bj. 
Following precisely the arguments given in Ref. 20 [see Eqs. (45)-(47)] we can then conclude 

N-l 



b =^~ 

3 p ' 



with 



i + E n# 



fc=i 



(71) 



(72) 



where the compact notation introduced in (8) and (15) has been used. 

To determine the coefficients dj and Tj that control the behavior of Cj(t) in (69) we substitute this ansatz into (68) 
concluding, for t — > oo, 



/•OO 

a 3 = { 1 (t) [ 0j -_i (t - r) + Tj_i] + p7 +1 (r) [a J+1 (t - r) + T j+1 ] 
Jo 

- K + W + ¥»7(r)] [aj(t - r) + Tj] - tpj +1 (r)b j+1 + ^(r)^} dr. 



(73) 
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Now we may introduce the first moments of the relaxation functions via 

roc 

- / T<pf(r)dT = (dy5f/ds)\ s=Q = gf, (74) 
Jo 

[see also (22)] which leads to 

a,j = t [uj_iaj_i + w j+1 a j+1 - (uj + Wj)aj] + [gf^aj-x + gj +1 a j+1 - {gf + gj)a 3 ] 

+ [uj-iTj-! + w J+1 T j+1 - { Uj + w 3 )T 3 } + [uj-xbj-! - w j+1 b j+1 ] . (75) 

The secular term proportional to t should vanish here, which condition requires 

= Uj-idj-i + Wj+idj+i — (uj + Wj)a,j, (76) 

while the coefficients Tj then satisfy 

a i = [ffZ-ifflj-i + 9j+i a j+i - (dt + 9j)aj] + [uj-iTj-i + w j+1 T j+1 - { Uj + w^Tj] 

+ [uj-ibj-i - w j+ ib j+1 ] . (77) 

Comparing (76) with (71) one can conclude that 

Oj = Ab 3l (78) 

where, using the normalization X^c/ = 1 following from (66) and (63), the constant A can be calculated by 
summing Eqs. (77) on j: the Tj cancel identically which leads to 

JV-l N-l 

Then, on using the result (72) for bj, we find 

A = N[l-n?]/R N . (80) 
To obtain the coefficients Tj we introduce, following Ref 20 [see Eqs. (54)-(57)], 

y j =w j+1 T j+1 -u j T j , (81) 

and rewrite (75) as 

Vj - Vj-i = a 3 - + 9j+i a 3+i - (.9/ + 9j)aj] - + w j+1 b j+1 . (82) 

The solution of this equation, which is achieved using the strategy described in Ref. 20, yields 

N-l 

Vj = u 3 b 3 + ( A / N ) + ^bj+i+i + i a i9t - a j+i9j+i) + c, (83) 

*=o 

where c is an arbitrary constant which will cancel in the final formula for the dispersion, D (sec Refs. 20 and 28). 
The fact that this expression solves (82) can be checked with the help of the relation 

Ujbj - w j+1 bj + i = A/N, (84) 

which follows from Eqs. (72) and (79). Then, iterating (81) and invoking the periodicity (70) yields the relation 



it,- 



N-l 

j+k 
3+1 

k=l 



(i-nf), (85) 



which, via (83), (78), (79) and (72), represents an explicit result in terms of the effective rates defined in (16). 

Now we can calculate the drift velocity, V, and the diffusion constant, D, using the long-time definitions (1) and 
(2). The mean position of a particle is given by 
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, +oo JV-1 , N-1 

<*(*)> ^EE^' + Ni)Pj(i,t) = - J2 cm 

l= — oo j=0 j=0 

With the aid of the generalized master equation (62) the derivative can be taken which, when t — > oo, leads to 

lim —(x(t)) = — {ui — Wj)bj = —A. 
t^oo dt /V v 3 3,3 N 



(86) 



3=0 



Using the result (80) yields our final formula for the drift velocity, namely, 



v a = d[i-n?]/R 



N- 



(87) 



(88) 



where we recall that Rn is defined in (15). This expression corresponds exactly to Dcrrida's original result for the 
simple sequential kinetic model. 

A similar approach suffices to determine the dispersion. We start from 



oo JV-1 



N 2 

1— — QO j — 

and again appeal to the master equation (62) in the long-time limit. This leads to 

N-1 N-1 



(89) 



3=0 3=0 



(90) 



Then, using (86), (87), and the definition (2), we obtain 
d 2 



D 



N 2 



N-1 N-1 N-1 N-1 

E K ~ w i) T i + \ E ("i + + E a M - 9j) -aJ2 t 3 

3=0 3=0 3=0 3=0 



(91) 



The coefficients Tj can be re-expressed using (85) and (83) from which the constant c then cancels 20 ' 28 . Finally, the 
definitions (19) and (20) allow us to write the dispersion in the form presented in (17)-(21) while the nonexponential 
parameters first introduced in (22) are confirmed by (74). Note, that the dispersion consists of two terms, D and 
D\, the second arising purely from the deviations of the waiting-time distribution functions from the "chemical," 
Poissonian forms. 



IV. PERIODIC MODEL WITH BRANCHES AND WAITING TIMES 

Now consider the one-dimensional periodic model with branches and waiting-time distributions as presented in 
Fig. 2. Let Pj,k(l, t) be the probability of finding the walker at site I in state j of the main sequence (labeled k = 0) or 
in state k = 1, • • • , L, on the associated side branch, at time t. Appealing again to Landman et al. 27 this probability 
is governed by the generalized master equation 

-p,, (M) = J {^tiMP^^oa^-^+^^p.+^oa^-^+^^^p.-^/^-r) 

'<pf(r) + <pJ(t) + ^ (r)] P jfi (l,t- r)} dr, (92) 



for k — 0, by 



dt' 



-P 



j,k 



(l,t) = {ip^^Pj^iht-^+^^P^+^t-T) 

'^)+'Pl k (r)\P j , k {l,t-T)}dT, 



(93) 



for 1 < k < L, while for k — L we have 
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tP, 



j,L 



(l,t)= f 

Jo 



dr. 



(94) 



dt 

The relaxation functions <pf(t), <fj(t) and <Pj(t), are defined, as before, via their Laplace transforms as specified in 
(6). At t = we may assume that Pj,k{h 0) = Pjk^lfi and normalization requires 

oo N-l L 



J = -00 j=0 fe = 



(95) 



Again following the strategy of Derrida's method as described in the previous section, we introduce the auxiliary 
functions 



B hk (t)= J2 p j,kQ,t), C jtk (t)= Y^{j + Nl)P^ k {l,t). 

I — — oo /— — oo 

The time evolution of £j,fc(i) is then described by the set of L + 1 equations 

d r* 

J t B i,o(t) = J o {^iWBj-i,o(i ~r) + <pJ +1 {T)B j+1 , {t - r) + p] A (T)B lA (l,t - r) 
- Vi» + VJW + ^.oW] s i,o(i - r)} dr, 
- ^ {^-iWBpt* ~r) + vlk + i{T)B ia {t - r) 
^, fc « + ^, fc (r)]B j ,o(t-r)}dr, 



(96) 



(97) 



(98) 



dr. 



fi,L-i<j)B itL _ x {t - r) - <pl L (T)B jiL {t r) 
Similarly, the time evolution of the Cj^(t) obeys the equations 

i,o(*) - / {^tiWCj-Lo^ - r) + ^7 +1 (r)C j+li o(i - r) + ^i(r)Cj,i(I,t- r) 
Jo 



+ ^"(r) + ^o(r)] Q, (t - t) + ^(r^-i^t - r) 



-Vj+i(T")-Bj+i,o(t-T)}dr, 
(*) = j[* - r) + ^ fe+1 (r)Q,2(t - r) 

-[^ fe ( T )+^ fe ( T )]c j -o(i-r)}dr, 



(99) 



(100) 



(101) 



d_ 

dt 



C jt L(t) 



i-iWCj,i-i(i - r) - <pl L (T)Cj, L {t r) 



dr. 



The previous arguments 20 ' 28 lead to the expectation 

Bj,k{t) — ► &j,fc, Cj,k(t) — a,j,kt — > 2},fe, 

for i — ► oo. At large times the equations of motion (97)-(99) then yield the relations 

Wj+ibj + i,o - Ujbjfl = wjbjfi - uj-ibj-ifi + {Pj,obj,o - Jj,ibj,i), 

Pj,obj,o - 7j,ify,i = /3j,i6j,i - lj,2bj,2 = ■ ■ ■ = Pj.L-ib^.L-i - 7j,Lbj,L = 0, 



(102) 



(103) 



(104) 
(105) 



where we have invoked the definitions (16) and (25) for the rates Uj, Wj, /3j t k and 7j.k- Recalling, likewise, the 
definitions (22) and (33), the coefficients and Tj.k must satisfy, first, 
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a 3,o = [gf-iaj-ifi + 9 j+1 aj+i,o - (.9/ + 9-j )aj, ] + [uj-iTj-1,0 + w j+1 T j+lfi - (uj + Wj)T jfl ] 

+ [uj-\bj-ifi - w j+1 b j+lfi ] + g] tl a jt i - gj fi a,jfl + 7j,i 7 j\i - Pj,oT jfi , (106) 
OjM = ko%,o + ffj, 2 a,,2 - (^i + .9J,iKm1 + [^,0^,0 + 7^,2 - (0 Jt i + 7^1)^,1] , (107) 



g^L^a^L-l ~ g],L a j,L + 0j,L-lT jtL -l - li,LT jtL 



while the vanishing of the secular terms yields, also, 

iOj+iOj+1,0 - MjOj.o = Wjfflj.o - Uj-iOj-i.o + Wj,oaj,o - 7j',i a j,i)) 

^j.offlj.o - 7j,i°j,i = ftM a j,i ^ 7j,20j,2 = • • • = 0j,L-ia,j,L-i - lj,La jtL = 0. 

The side-branch functions frj^ are found easily by solving Eqs. (105) recursively which gives 



0j,k 



nf%, (fc = i,---,L), 



(108) 



(109) 
(110) 



(111) 



where we have used the product notation (9). To obtain an expression for bjp we follow exactly the method used to 
derive (72) in Section III thereby finding 



bj,o = rj/R^, 



(112) 



where rj is defined in (15) while Fi N was introduced in (24). Comparing (104) and (105) with (109) and (110) leads 
to 



a jik = Ab jik with A = N(l - U^)/R P N . 



(113) 



The final derivation of the expressions for the drift velocity and the dispersion now follows along the lines developed 
in the previous section for the models without branches. The results have been presented in full in Section LB. As 
regards the velocity the branches generate no changes beyond the replacement of rj by r^ and Rn by R^ N . However, 
an additional term, D 2 ^ 1 appears in the dispersion: see (31)-(33). 



V. PERIODIC MODEL WITH DEATHS AND WAITING TIMES 

Consider, finally, the periodic sequential model with waiting-time distributions and the possibility of an irreversible 
detachment or death from each state that is described by a waiting-time distribution function ipj(t) — ipj ±N (t): see 
Fig. 1. The generalized master equation for the probability Pj(l,t) now reads 27 

j t P 3 {l,t)=J { <Pt-i WPj-i (l, t-r) + (r)P j+1 (I, t - t) 

[<p+(r) + <fij(r) + ^(r)] Pj{l,t- r)} dr, (114) 

where, as before, the relaxation functions <pf(t) and <Pj(t) are related to waiting-time distribution functions via (5) 
and (6). We may again assume that the initial condition is 



N-l 



Pj(l,0) = Pj\ with P J = L 

3=0 



(115) 



However, as discussed in Ref. 20, because the total probability is no longer conserved [so that Y^t^oo S^Lo 1 ^'(^ * > 
0) < 1], we look for long-time solutions of the generalized master equation (114) that are of the form 

P,(M)«e- At -^(M), (116) 

where, as before 20 , the decrement and the periodic state coefficients, tj = tj±n 7 are to be found from the requirement 
that Pj(l,t) satisfies a suitably "renormalized," probability conserving master equation 
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j t m,t) = u j _iP J -_i(i,t)+s J -+iP j+ i(j,i) - +wj)Pj{i,t). (117) 

By substituting the ansatz (116) into the full generalized master equation (114) we obtain, for large times, the equation 
j t Pj(l, t) = Uj-^-^Pj-! + w j+1 e T *-^P j+1 - (uj + Wj + 5j - X)Pj, (118) 

in which the modified rate definitions (35) et seq., which depend explicitly on A, have been used. Matching terms 
with those in (117) generates the identifications 

Uj—UjEj+x/sj and wj = UjEj-i/sj with ej — e Tj . (119) 

It also yields a condition which the e 3 - must satisfy for consistency, namely, 

- WjEj-i + [uj + Wj + 6j)ej - UjS j+1 = Xej, (120) 

But, recognizing the periodicity in j, this is precisely equivalent to the eigenvalue equation Me = As, where M is 
the NxN matrix specified in (12) and (13). Since the asymptotic decay is required in (116), A must be the smallest 
eigenvalue which, clearly, should be real and positive. 

To find expressions for the drift velocity and the dispersion we now require three auxiliary functions, namely, 



oo 



Bj(t)= J2 W), Cj(t)= £ U + Nl)W,t), (121) 

I— — oo l= — oo 

and also 

oo 

Ej(t)= E U + NlfPj(l,t). (122) 

l— — oo 

For large t we may expect the asymptotic behavior 

Bj(t) ta e- x *- T > Bj(t), Cj(t) « e" A *-^ Cj(t), and Ej(t) » e- A *- T '25,-(t), (123) 

with, extending Derrida's ansatz, 

Bj(t)^bj, CjW-ajt^Tj, and E^t) - ejt 2 - fjt - X s . (124) 

The explicit formulas for the coefficients bj , cij , Tj , ej , fj and Xj can now be found straightforwardly by extending 
the procedures outlined in Sections III and IV. However, the detailed calculations are fairly tedious and, because of 
the presence of the functions Ej(t), give rise to the higher order nonexponential parameters, and hp defined in 
(47) and (48). 

The mean displacement at time t must now be suitably normalized so as to include only surviving walkers. Thus 
we have 

, oo JV-1 / oo N-l , N-l I N-l 

^W) = n £ E(H^(M) /EE p j( 1 ^ = n E cm / E B M (125) 

J=-oo j=0 I /=-oo j=0 j=o I j=0 

while the desired mean-square displacement is similarly given by 

, 2 oc JV-1 / oo JV-1 , 2 N ~ 1 /JV-1 

^ 2 (*)) = N2 E Ec?+ M ) 2 ^(M) /EE = ^2 E ^-(*) / E ^w- ( 126 ) 

J=-oo j=0 / ;=-oo j=0 j=0 I 3=0 

On taking the derivatives required by (1) and (2) and the steady-state limit, these expressions yield the results for 
the velocity and the dispersion given in (34) and (39)- (48). Naturally, when there is no possibility of detachments 
[ipj(t) = 0, A = 0] one recovers all the results for the periodic sequential models with waiting times as reported in the 
Introduction and Summary and in Section III. 

This completes the description of our mathematical analysis. In brief summary, we have introduced linear, peri- 
odic sequential stochastic models with general waiting-time distributions and have found explicit expressions for the 
corresponding mean velocities and dispersions that have been reported in Section I. The simplest sequential models 
have been extended by including finite branches and by allowing for the possibility of death or detachment processes. 
The deviations from the exponential waiting-time distribution functions that characterize standard kinetic models 
embodying Poisson processes, do not change the form of the velocity expressions; however, the dispersions entail 
nonexponential parameters that enter in a more complicated manner. The concept of "mechanicity," introduced in 
Section II, is useful to quantify and visualize the departures from the usual "chemical" kinetic descriptions. 
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TABLE I. Expressions for the Rate Moments in Terms of Waiting-Time Moments. See Eqs. (5)-(7) and note 
that vi = for n > 1 when ip c {t) oc ip(t) oc e~ ct for any c (> 0). 
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Fig. 1 Kolomeisky and Fisher 



FIG. 1. A schematic general linear, periodic, stochastic process with forward and backwards waiting-time distributions, 
ipj~(t) and 4>~{t), failing or futile attempt distributions 4>j(i), and irreversible death rates, ipj(t), from states j t , where the 
reference states, Oj, are located at positions x = Id. (Note that the precise locations of the intermediate states 1;, 2;, • • • , (N — 1);, 
have no significance for the velocity, V, and dispersion, D, since these are defined asymptotically for large times, t — ► oo.) 




Fig. 2 Kolomeisky and Fisher 



FIG. 2. A schematic periodic stochastic process with branches of finite length (< L) grown from each primary site 
(j, 0)( = ji; outward and inward waiting-time distributions, ipj k (t) with k = 0, 1, • • • , (L— 1), and ipl k (t) with k — 1, 2, • • • , L, 
are specified at each branch site (j,k)i. Failing or futile attempt distributions, %f)° k {t), are not shown but may be considered 
as present. 
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